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Abstract 
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1 Introduction and Summary 



One of the most interesting outcomes of the AdS/CFT correspondence [1, 2, 3, 4] is the 
abihty to study quantitatively the deconfined phase of 3+1 dimensional gauge theories, 
something which cannot be done analytically for QCD (except at temperatures much 
higher than the QCD scale). For strongly coupled large N gauge theories which have 
(at zero temperature) a dual description given by a weakly curved string background, 
the deconfined phase has a dual description in terms of a black hole (black brane) back- 
ground which can be reliably studied in the supergravity approximation. Of course, 
theories with a weakly curved dual are rather different from QCD in various ways, but 
one can still hope that their deconfined phase will not behave that differently from that 
of QCD, and in some cases these theories are even continuously connected to (large 
A^) QCD by varying a dimensionless parameter, and one could hope that the depen- 
dence on this parameter is small (at least at temperatures of order the deconfinement 
temperature). 

The simplest theory to study in this way, on which most of the research thus far 
has focused, is the strong coupling limit of A/" = 4 supcrsymmctric Yang-Mills theory, 
whose deconfined phase has a very simple description as a black hole in anti-dc Sitter 
space [5]. This theory does not confine at low temperatures, but its deconfined phase 
still seems to exhibit many similarities to that of QCD. Obviously, it would be nice to 
have additional examples of 3 + 1 dimensional deconfined theories which can be studied 
quantitatively by using their gravity dual, and, in particular, examples of deconfined 
phases of confining theories, in which one could study the dependence on the temper- 
ature compared to the deconfinement scale^. So far there is only one known example 
of such a deconfined background, which is that corresponding to 4 + 1 dimensional 
supersymmetric Yang-Mills theory (with a specific UV completion) compactified on 
a circle with anti-periodic boundary conditions for fermions [5, 6]. It would be nice 
to have additional examples, especially since in the example above the physics at the 
deconfinement scale is really five dimensional rather than four dimensional. 

In this paper we study the deconfined phase of the confining "cascading gauge the- 
ories" constructed by Klebanov and collaborators [7, 8, 9]. The equations determining 

^Notc that in any large N gnugc theory with a weakly curved string theory dual the deconfinement 
transition is a first order phase transition. This is similar to large N QCD, but it is different from 
QCD itself, so one should be careful when comparing the behavior of such theories near the phase 
transition to that of QCD. 
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the corresponding black hole solutions are quite complicated, and have no known an- 
alytic solutions. At very high temperatures it is possible to find analytic solutions in 
an expansion in inverse powers of the logarithm of the temperature, and the leading 
order solution in this expansion was found in [10]. This solution shows that at high 
temperatures the number of degrees of freedom in the theory grows as the square of the 
logarithm of the temperature [11, 12, 10]. In this paper we numerically solve the equa- 
tions for a wide range of values of the temperature, in the supergravity approximation, 
and use the solutions to analyze the thermodynamics of the deconfined phasc^. 

At low temperatures the "cascading gauge theories" spontaneously break a dis- 
crete chiral symmetry (and also a continuous U{1)b symmetry) [9, 14, 15, 16]. At 
high temperatures one expects these symmetries to be restored, but apriori it is not 
obvious whether there is a single phase transition from a low-temperature phase with 
confinement and chiral symmetry breaking to a high-temperature phase with no con- 
finement and chiral symmetry restored, or whether additional phases also exist. In this 
paper we only study deconfined phases in which the chiral symmetry is restored; we 
plan to study the possibility of having other phases in the future [17]. In the classi- 
cal supergravity approximation the free energy of the low temperature confined phase 
vanishes (since it only arises at one- loop), and thus the phase transition to a black 
hole background occurs at the lowest temperature for which the free energy of a black 
hole background starts becoming negative. Assuming that this transition goes directly 
to the chirally symmetric black holes that we construct, we find that there is a first 
order deconfinement transition at a temperature TcnUcai = 0.614111(3)A, according to 
a specific definition of the strong coupling scale A that we describe in section 5^. The 
black hole backgrounds continue to exist also at lower temperatures, but they have 
positive free energies so they no longer dominate the thermodynamics. Presumably, 
as the temperature is lowered further, the black hole backgrounds eventually become 
singular; in this paper we only compute the numerical solutions until a temperature 

^Similar numerical solutions were studied in [13], but we do not understand the parametrization 
used there to analyze the solutions. Presumably, our solutions should be identical to (some of) the 
solutions of [13], but our parametrization allows for a direct computation of the thermodynamical 

properties of the solutions. 

^It is easy to translate this definition to other definitions of the strong coupling scale, such as the 
mass gap. Note that, as in all theories with a gravity approximation, the square root of the confining 
string tension is not a useful measure of the strong coupling scale, since it must be much larger than 
all other measures of this scale for gravity to be a good approximation. 
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slightly below the deconfinement temperature, so we do not see this. 

There are several interesting directions for further study. We are currently working 
on checking whether the solutions we find are stable to deformations which break 
the chiral symmetry, in order to sec if there are signs of a deconfined non-chirally- 
symmetric phase appearing at intermediate temperatures [17]. The black hole solutions 
that we find (numerically) can be used for a detailed analysis of the properties of the 
deconfined phase (for instance its hydrodynamical properties [18] or jet quenching 
[19, 20]); it would be interesting to understand the similarities and differences between 
these properties and those of deconfined conformal theories. It is possible to add fiavors 
in various ways to the "cascading gauge theories" [21, 22, 23, 24, 25, 26, 27], and it 
would be interesting to study the flavor physics in the deconfined phase, and whether 
there are any phase transitions associated with the flavor sector. 

Our study is purely in the supergravity approximation; it would be interesting to 
study the corrections to this approximation, in particular those coming from string 
theory corrections to the supergravity action. The cascading gauge theories have a 
continuous dimensionless parameter such that in one extreme of this parameter super- 
gravity is a good approximation, while in the other extreme they reduce to a standard 
M = 1 supersymmctric Yang-Mills (SYM) theory. At zero temperature supersymmetry 
tells us that the dependence on this parameter is smooth. However, it is not obvious if 
the behavior of the deconfined phase is smooth as this parameter in changed; in fact, it 
seems plausible [28] that as in other similar cases [6] there would be a phase transition 
in this phase, since the geometry in the supergravity regime does not have any cycles 
shrinking at the horizon (except for the thermal S^), while in the SYM regime one 
expects the transition to be independent of the KK modes so a two-cycle should still 
shrink (as it does in the confined phase). Of course, even in the absence of such a phase 
transition, the behavior of the theory could be modified as the dimensionless parameter 
is changed, so our analysis does not teach us directly about the phase structure of the 
large N U = I SYM theory. 

This paper is organized as follows. In section 2 we describe our ansatz for the black 
hole solutions and the equations of motion that it leads to. In section 3 we describe 
the boundary conditions for these equations. In section 4 we analyze the meaning of 
the parameters appearing in these boundary conditions, and show how to translate 
them into physical quantities such as the temperature and free energy. In section 5 we 
describe our numerical procedure and present the "bare" numerical results. In section 



4 



6 we translate these results into physical properties, and present the results for physical 
quantities such as the free energy and the expectation values of various operators as a 
function of the temperature. An appendix contains a pcrturbative analysis of the very 
high temperature solutions; this is useful both in order to make sure that our analysis 
is valid by verifying that it is consistent (at least at very high temperatures) with the 
first law of thermodynamics, and in order to test our numerical solutions at very high 
temperatures by testing their agreement with the perturbative expansion. 

2 The equations for the cascading black hole 

In this paper we compute the metrics corresponding to the finite temperature behavior 
of the "cascading gauge theory" found in [7, 8, 9]^, which may be thought of as a specific 
SU{K) X SU {K + M) M = 1 supersymmetric gauge theory, with a number of colors 
K which runs logarithmically with the energy scale [9, 11, 12, 30, 31]. The "cascading 
gauge theory" has a single dimcnsionless parameter (in addition to the integer M). 
which in the gravitational description of this theory can be taken to be QsM where 
Qs is the string coupling (which is constant in the zero temperature solution) and M 
is the RR 3- form fiux (corresponding to the number of fractional brancs). When this 
parameter is large, the gravitational description of the background is valid at all scales. 
On the other hand, when it is small, the theory at low energies reduces exactly to the 
M = 1 SYM theory, but the gravity dual is highly curved. We will only analyze the 
theory in the regime of large QsM ^ where the gravitational approximation is good and 
all radii of curvature are large compared to the string scale. 

As in any other confining background, the low temperature behavior of this theory 
is governed by a gas of hadrons; the gravity dual of this description is simply given by a 
thermal identification (t = t+l/T) of the zero temperature solution found in [9]. As the 
temperature is increased one expects the theory to deconfine; in the gravitational dual, 
deconfined phases are described by black holes (whose horizon fills all of space, so they 
are really black branes). Our goal in this paper will be to compute the gravitational 
backgrounds corresponding to the deconfined phase of the cascading gauge theory. Note 
that the low temperature phase is stable all the way up to the Hagedorn temperature 
of the confining theory (related to the confining string tension); when the gravitational 
approximation is valid, this temperature is very large compared to the characteristic 

*See [29] for a recent review of this theory. 
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mass scale of the gravitational background (which determines the mass of the low-lying 
hadrons). Thus, as in all other cases of confining backgrounds with gravity duals, we 
expect the deconfinement transition to occur at a temperature which is much smaller 
than the Hagedorn temperature, which means that it should be a first order phase 
transition. 

The cascading gauge theory has a Z2M chiral symmetry [9, 15] which is sponta- 
neously broken to Z2 at low temperatures (by gaugino condensation in the limit where 
the theory is a pure SYM theory), and it has a U{1)b symmetry which is also sponta- 
neously broken [14, 16]. At high enough temperatures we expect these symmetries to 
be restored [11, 32]; this expectation is confirmed by the analysis of the asymptotically 
high temperature black hole solutions in [10]. Apriori it is not obvious if the decon- 
finement transition happens together with the global symmetry restoration transitions, 
or if the transitions are separate. In this paper we will only look for solutions which 
preserve the chiral symmetry and the ^7(1)^ symmetry; the stability of these solutions 
with respect to chiral-symmetry-breaking deformations will be analyzed in [17]. We 
also assume that the solutions preserve the SU{2) x SU{2) global symmetry of the 
theory, which is preserved also at low temperatures so it is reasonable to assume that 
it is preserved at all temperatures. 

The form of the gravitational background of the cascading gauge theory at large 
radial variables (close to the boundary) was found in [8] and is known as the Klebanov- 
Tseytlin solution; this form preserves the full global symmetry. The solution at any 
temperature is expected to asymptote to this background near the boundary (the "UV 
region"). In the zero temperature solution [9] the Z2M and U{1)b symmetries are 
broken far from the boundary (in the "IR region"), but we will look for solutions 
where they are preserved. Recall that the solution for M = (where the theory 
does not cascade) is [33] AdS^ x T^'^, and that the solution of [8] has a similar form 
but with the radii of curvature (and the flux) varying logarithmically in the radial 
coordinate. Since we are looking for solutions that preserve the full global symmetry, 
we can perform a Kaluza-Klein reduction on the T^'^, and leave only the flelds which 
are singlets of the global symmetry group. In general there are 5 such fields^; the five 
dimensional graviton and 4 scalar fields. In the M — Q limit the scalar fields have 
scaling dimensions A = 4.4.6.8. The scalar fields are various linear combinations of 

^We only consider the fields which are turned on in the solutions that we are interested in; other 
fields, such as the type IIB axion, are consistently set to zero. 
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the dilaton, the overall volume of the T^'^, the relative size of the circle in T^'^ (thinking 
of T^'^ as a circle fibration over two 2-spheres), and one mode coming from the RR 
fields. 

We are looking for black hole solutions that preserve spatial rotational and trans- 
lational invariance, as well as time translation invariance, so we can always choose a 
form of the five dimensional metric where only 3 components are not vanishing - Gqq, 
Gii {i — 1,2,3) and Grr (where r is the radial position). We can use the freedom of 
reparametrizing the radial coordinate to eliminate one of these degrees of freedom - 
we will choose our radial coordinate x to be defined by 

p=-{i-xr (2.1) 

(with no summation over i). This choice is convenient since at the boundary we expect 
the metric to be Lorentz-invariant so a: — > 0, while at the horizon Gqo vanishes so 
a; — > 1; the range of the radial coordinate in our parametrization is thus x G (0, 1). 
Note that in the conformal M = case there is a simple form of the black hole solution 
in this parametrization, which is given by 

dsl, = e (2X - X^)-V2 _ ^)2^^2 ^ ^^2 ^ ^^2 ^ ^^2^ ^ dx^ ^ (^^1,1)2 ^ 

(2.2) 

where (rfT^'^)^ is the metric on T^'^, and the constant ^ is related to the temperature 
T as follows: 

e = vrT. (2.3) 

Motivated by the form of (2.2), we write down the most general ansatz for a black 
hole metric preserving all the symmetries as^ : 

a=l 

where h, /2 and arc some functions of the radial coordinate x. There is also a dilaton 
g{x), and form fields given by 

K 

Fg = P ey, A (e^i A e^^ - ee^ A e^J , -B2 = — (e^^ A e^^ - e^^ A e^J , 

(2.5) 

-P5 = -^5 + ^-^5 , = -K A A e^,^ A ee^ A e^,^ , 



^Thc frames {eg^, e^^} are defined as in [30], such that the metric on a unit size T^'^ is given by 
(4)+ELi K+e^J. 



where K is a function of the radial coordinate x. The constant P appearing in (2.5) 
is a constant times the quantized flux M; we will write down the precise constant in 
terms of the five dimensional Newton's constant below. We will find it simpler to work 
in terms of P rather than M, and we will only go back to using the integer flux M 
in the final section. After we gauge-fixed the radial coordinate by (2.1), we have a 
constraint equation coming from the equation of motion of this variable; we can use 
this equation to solve for Gxx, which is given by 

= 2(.-l)P$g-.)^.^A (l2^VI//./.^(l - .) 

+ hx\2P'f,g'h" - UPYh'f!,' + K"hg + 2P'h'f,g")ix - 1)(2 - xf 
- AxP'hg'hM2 - x)ix' -2x + 2){h' h + 

+ 4xPV3(2 - x)g-'h\2xf^{l - x){2 - x) - [x" -2x + 2)U)f^ , 

with 

A = K' + Sh^flh{f2 - 6/3) + 2hflP^g. (2.7) 

All in all, we have 5 scalar functions of x that we need to solve for : h, f2, fs, g 
and K. We can derive the equations of motion for these fields, in the supergravity 
approximation, either directly from the ten dimensional type IIB supergravity action, 
or by first reducing this action to five dimensions and then deriving the equations of 
motion. The equations that we find take the following rather complicated form : 



0^h"-[8hf2{f2-6f3) + gP' 



' A 

2 d2. 



+ [8x/iV3V2(/2-6/3)(x-2) 
h' 



- 6{K' + hfiP'g)^ + [%h'flh{h - 6/3) + + AhflP'g] 



K 
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+ K\^x'-Qx + A)+AhftP'g{l-xY]-— 

x{l — x){2 — x)I\ 

'/|A 

+ {K' + hflP'g)^ - 2h{K' + hflP'g)[2xai - x){2 - x) - {x' - 2x + 2)/3] 



X 



/3/2x(l-x)(2-x)A 



+ 8(x' -2x + 2){K' + hfiP'g) 



4gfiP^A 
[X-- 



x(l-x)(2-x)fsA 



2[7hfiP'g + 16h'fif2{f2 - 6/3) + 5K' 



h 



{2-xyx^A' 



(2.8) 
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- ^ - [4x/^/3/i(l -x){2- x){gP' - 8/|/i) + Shh'fUdxMx - 2) 

J2 



f 

+ /2(4 - 2x + x^)} + a;i^2(2 - a;) - ihfiP^g{l - x)^' ''^ 



{1 - x){2 - x)xA 

hA 



fl{gP' - 8f,h)^ - 6h{gP' - Sf.h)^ 



- {K^ + 2Afih%ih - m) Jfp,g^ + hfiigP' - (2.9) 
+ 2/32(^2 -2x + 2){gP^ - 8fih)- 



{2 - x){l - x)xA 
+ 8hf3{x'-2x + 2){gP'-8fih) 



52 0/27N 72/3 



{2 - x){l - x)xA 
+ 2[16fih'M2h - 3/3) +K'- hfiP'g]- 



(x- 2)2x2 A' 



=/^' - [2P'ghfi + 8flhh\Ah - 15/3) + ^ih - S/s)^^^ 

+ 2(/2 - 3/3)^^^ + Kh - 3/3)-^^'^'^^'" 



P2^A ^2 A 

/^/|/2 

(2-a;)(l -x)xA' 



+ 8(x^ - 2x + 2){h - miT. . , h' 



- 8{h - 3/3)M(l -x){2-x)- {x' -2x + 2)/3] (^3^y^^#^ ^^'^^^ 
+ {8/2/|[(5a;' - lOx + 8)/2 - 6/3(3a;2 - 6x + A)]h'' + 2xP^gfl{x - 2)h 

+ "^^'^"^ " ^)> (1_^)4_^)A + + - ^^^'^^^^^^ + 

X 



x2(2-x)2A' 



KK'"^ K' 
=ii'" - — ^ + N/3 + 2/i(?/^(l -x) + hhg'H -x)+ ghh\l - x)]- 



^ ,2P2^ _ _ 2P2^ + 4[2./^(l - .)(2 - X) 

- (x^ - 2x + 2)/3] ^ , KahP^hf^ _ 4p2^^(^2 _ 2^ + 2). ^3^' 



72(l-a;)(2-a;)a;A ' ' (2 - x)(l - x)a;A 



- 16(x2 -2x + 2)— ^ + 12- 



(2-a;)(l -a;)a;A a;2(2-a;)2A' 

(2.11) 
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- ^ - [8/iV3V2(/2 - 6/3) + W'ghfl + K']^- P'^-^ + 6P'^ 
+ {^h'flHh - 6/3) + i^^)T^Srx - 2/l(^^ - 2x + 2)P- 3'^' 



2hiP\2xr^{\ -x){2-x)- {x' -2x + 2)h' 



f2{2-x){l-x)xA 



8(x 2x + 2)P (2-x)(l-x)xA+^^ x2(2-x)2A- 

(2.12) 

3 Boundary conditions 

In order to solve the equations of motion (2.8)-(2.12) we need to specify boundary 
conditions, both at the asymptotic boundary and at the horizon. We will require that 
asymptotically the solution should match onto the Klebanov-Tseythn (KT) solution, 
and that it should be regular near the horizon. 

3.1 The UV boundciry conditions 

Near the boundary a; — it is possible to solve the equations by a power series in x 
and ln(a;), whose leading term gives the KT solution. This expansion takes the general 
form : 



2 00 n 

h =Vo - ^ H^) + E E ^--'^ 

oan — — 



n=l k=l 



/2 =«0 + E E ^""''^ ^"^^ 
n=l k=l 

00 n 

h =«o + E E ^-.^ ^^'(^)' (3-1) 

n=l k=l 

00 n 

K ^Aho,^al - -P'go - -P'go ln(x) + E E ^-,k x^'^ In^(x), 

n=l fc=l 

00 n 

^ + E E ^"''^ 



n=l fe=l 



Most of the coefficients appearing in this expansion are not independent; the indepen- 
dent coefficients correspond either to parameters of the cascading gauge theory or to 
vacuum expectation values (VEVs) of the operators dual to the fields we are solving 
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for. In the KT case there are 3 asymptotic parameters, which we choose to be go, 
/io,o ^-nd ttQ. go is related to the dimensionless parameter of the cascading gauge the- 
ory, and we will see that one combination of the other parameters is related to the 
temperature and the other is related to the dynamical scale of the cascading theory. 
Naively we would expect to have 5 parameters related to VEVs, but in fact there is 
one relation between the VEVs which is given by the conformal anomaly equation (see 
[30]), so we are left with four parameters corresponding to VEVs, which we choose to 
be {a2fi, 92,0,(13,0, (i4:fl}- Note that a VEV appearing at order x"/^ corresponds to an 
operator which has dimension 2n in the conformal limit of the theory. 

Using these 7 parameters we can solve for the coefficients in (3.1) to any order 
we wish. It turns out that there are no non-zero coefficients at order 0(x^/^). The 
non-zero coefficients at the following orders are : 

■ order 0{x): 

, _ ^90^2,0 p2 h - ( ^900-2,0 90 92,0 \ p2 , 3a2,o/io,o , . 

^^■^-"^8^^' 1,^8^" IK" IK J 

, 1 65(002,0 „s 

02,0 = -^^02,0, 5^2,1 = rj^^ , (3.3) 

6^002,0 p2 ( 90^2,0 1 1 \ p2 , 24 

K2,i = P , K2,o = [—^ -go - 2^2,0 ) ^ + y«oa2,o/io,o- (3.4) 

■ order 0{x^''^): 

^^■°-"60^^' ^3,o--^a3,o, ^3,0--^ P. (3.5) 



Order 0{x'^): using the notation 6 = 139P^go — 120/io,oflo' '^^ have 



,2 io„2 io^2 



'^'^ ~ 196a^ ' ' 245an ' 245a, 



3^002,0 „2 _ 12ai o , _ 12a 



h.,3 = -"-^P\ «4,2 = -^, ^4,2 = -^, (3.6) 



"0 

, ^1[/ 142637gga|o 75^^02,0 S^g 75^^04,0 5c/|o 139goQ2,o^2,o ' 
^'^ 5\V 3920a^ 16a3 16al Sa^ 16a§ 280^ 

1002ffo/io,oa|o _ 1440/igoa|o 
49ag 49 



(3.7) 
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+ 
+ 



95gQ 230383gfoa2^o _ 139^oa2,o5'2,o _ 35goa4,o _ ^92,0 _ 1423goa2,o \ 4 



32a2 



23520a^ 



24al 



224a3 



/ 1363275(o/io,oai,o , 5 5/io,ofi'|o 528/io,oa2,o5'2,o 1505foa4,o/io,o 

o^r 2 ^ nQoii-ofi 

245a5 2 g^o 7ao ao 



570go/i,o,oa2,c 
7ao 



5160/ig,oa|o 
49 



+ 



5760alhloalo 
49^0^' 



(3.8) 



1 I ^ 31973^^a^^o ll^ga4,o 201c/ga2,o _ 219c/g _ 139^o^2,o 695c/oa2,oc/2,o 

~r <q ~r i-^q ^^^o ~r 



17640a^ 



56al 



32al 



252al 



288ag 



f 355foa4,o/io,o , 675fo/io,o«2,o 1675fo^o,o 262231c/o/io,oai,o 



+ 



V 2ao 
28/io,o02,ofl'2,o 15 



+ 



4ao 



2940ag 



5/io,o5'|c 



+ -rfl'2,0^0,0 + „ 
4 3t/o 



+ 600aoa4,o/io 



^ 2104ao/io,o«2,o^2,o _ 535356/ig_o«2,o _ 2280aoa2,o/io,o ^ 20ag^|o^o,o 



^90 

49536ag/ij],oa2,o 



245 



^0^ 



04,1 



49^0^2 

12463yoai,o 15 



70a, 







2 

1 '^ofi'2 

- -^5'o«2,o - ttS'oOo + 155-004,0 + ^ ' 
z z zgiQ 



26688 2 , 48a2/io,oa2,o52,o USSa^/i^ ^ai^ 

H O'0i^0,0O'2,0 "I ^ 



245 



7^, 







49P2^o 



1 

+ 



7177^0^2,0 15 1 , , ao5i,o 

— TTTT^ -?r9oa2,o - 77^0^0 + 15c/oa4,o + -7, 

490ao 2 2 2gQ 



480^/10,002,0^2,0 1152alhl 



0"2,0 



75'o 



49P2^o 



ir/6366^ 74 

^1 i 245ao ^ ° y5'o«2,o - 415-0^4,0 



(3.9) 

p2 

(3.10) 

7872ao/io,o02,o 
245 

(3.11) 



p2 ^ 126144ao/io,oai,o 



p2_ 



576ao/io,oa2,o5'2,o 
75'o 



^0 y 245 

,31,2 „2 



120aoa4,o/io,o + 



480a2,o/io,oag , 13824a^/i4o«2,o 



+ 



49P2^o 



9^,2 



ISfi'oo.i^o 
49a2 ' 



54,1 



365002,0 , 35002,0 , 602,052 



49a2 



+ 



7an 



+ 



7an 



(3.12) 
(3.13) 
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1 r /214407(?o'a^_o , , 32435o'«2,o 225(/2a4,o , 263^0 , 139^2,0^70^ 



^^'""^tV 980a2 + 4 ^ 28ao 2ao ^ 4 ^ 

7200^o^o,oai,o 360ao/io ,o«2,ofi'2,o 360ao5'o/io,o«2,o 2 7 

H |g y DUao52,o/io,o 

eOag^o^o.o ^ 8640ag/ig_o«2,o 



50 49P2 



(3.14) 

^--^^^ (3.15) 



"0 

Iff migla2,o 9719.r/2«2o 195a2a4,o , l^lo ISg^ 



5 W 56ao 56aQ 4ao 8 8 

417g'oa2,og'2,o Ap4 / 360fl'oaoa2,o/io,o 51096g'o/io,oai p 

14ao /V 7 245 ^ ' ^ 

, , 12384aLa2/i2, 

+ 48ao/io,oa2,o5'2,o P 
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189^^ 123^^02,0 , 97^^04,0 , 23885(?o'a2,o , 139^002 ,o5'2,o 139g'o5'2,o 

2 




^'°~(^\V 16 56ao 4ao 392a2 ^ 14a, 



57^1,0 ^ p4 ^ |^45c/oaoVo , 283077c/o/io,oa2,o , on„2 







, - . , ^ H ^ttH — ^ + 30005^2,0^0,0 - 3905oOoa4,o/io. 

4 y \ 1 z4o 

_^ 29735oaoa2,o/i-o,o _^ 240ao/io^oa2,off2,o \p2 _ 1440a[j/ig 002,0^2,0 



(3.17) 



75-0 

1440a^/io o«2,o 3571202 oOo^o 



+ 
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3.2 The IR boundary conditions 

Next, wc discuss the behavior of solutions to (2.8)-(2.12) near the horizon, a; — >■ 1. 
Introducing a near-horizon coordinate 

y=l-x (3.18) 

we find that in order for the solutions (2.4) to have a non-singular Schwarzschild hori- 
zon, the functions {/i, /2, /s, K} must all be even functions of y with a good Taylor 
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series expansion around y — : 

oo oo oo 



2n 



n=0 n=0 n=0 

oo oo 



(3.19) 



n=0 n=0 



When one solves the equations of motion perturbatively in y, one finds that the 
solutions are labeled by six independent parameters, which one can choose to be 
{Hq, Oq, 6q, /cq , Qq, Qi}. Naively one might think that all five equations of motion (2.8)- 
(2.12) would have one normalizable mode and one non-normalizable mode near the hori- 
zon, so that requiring a regular solution will set the coefficients of the non-normalizable 
modes to zero and leave us with 5 parameters. However, it turns out that for one com- 
bination of the equations both modes are normalizable near the horizon, leading to the 
extra parameter; this is related to a scaling symmetry of the geometry (2.4) which we 
will discuss in the next subsection, which implies that one combination of parameters 
is not determined before choosing a scale ^. 

Using these 6 parameters we can solve for the coefficients in (3.19) to any order we 
wish. Using the notation 5'^ = 8/io(flo)^ ~ P^9o: we have 
■ order 0{y'^): 



&t = ^ \)^gXP' - 2/^0^0 (3(aS)' - KK + K'A - ^0.%) } , (3.20) 



'^i - Aft 1 r ' ^1 



(5M j ' 2ag 



order Oiy'^): 

= Wf {- ''°'"f°° + "'»' (4ftS(aS)^ - UIJM - siP') ] . (3.21) 
«5 = {- ''"''"^ilj)^'"''' - <4yi - siP') } . (3.22) 



''Such a scaling symmetry in geometries with translationaUy invariant horizons was also noticed in 
[34]. 



14 



(3.23) 



(3.24) 



((5'^)2 

+ 180(6^)^4a? - %AaX{a\f + 180(6^)^(a?)^)6o^ + (6(ao^)^ + 9foo^(ao^)^ 

+ I2ala% + 8a^(a^)^ - 12(a?)^6^) - \{g^)%P''^ . 

(3.25) 

4 Mapping of parameters to the field theory 
4.1 Translation to the pcirametrization of [30] 

In this section we wish to understand the physical meaning of the 3 parameters which 
we used in the previous section to parameterize our theory (in the UV) - /io,0) ^-nd Qq. 
As we mentioned, g^P is the dimensionless parameter of the cascading theory (which 
must be large for the gravity approximation to be valid), while /iq.o and are related 
to the scale of cascading theory and to the temperature. Note that our ansatz (2.4) is 
invariant under a scaling symmetry taking 

{t,x)^\-\t,x), h^X-'h, /2^A/2, /3^A/3, (4.1) 
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and leaving all other functions in our solution (as well as the coordinate x) invariant. In 
terms of our asymptotic parameters, this scaling transformation leaves /io.oc^o invariant, 
meaning that this combination is a function of the dimensionless parameter of our 
theory, which is the ratio between the temperature and some scale A which characterizes 
the cascading theory. We can choose this scale A to be, say, the mass of the lightest 
glueball, or the square root of the string tension. We will find it more convenient to 
use a different definition of A which will be described below. 

The parametrization (2.4) we used above for the solution breaks down at zero 
temperature, since it assumes the existence of a horizon. In order to understand which 
combinations of our parameters depend on the temperature and which do not it is 
convenient to switch to a different parametrization of the geometry, which is valid also 
at low temperatures : an example of such a parametrization is given by [30] 

ds?„ =A-"V-' (-/"<<*" + dxl + dxl + dxt) + h}"'p-\d(,f 

a=l 

where {h, /2, /s, /} are functions of p. In this parametrization the supersymmetric 
zero temperature solution is characterized by two parameters: the value of the string 
coupling go, and the coefficient of the warp factor ho,o; in terms of these parameters 
we can write the asymptotic solution for small p as 

h = ho,o - hoP'' Hp) , K = 4/io,o - hoP^ - 2^o^' Hp) , 

2 ^ (4.3) 

g = go, / = A = /a = 1 • 

Note that the ansatz (4.2) is invariant under a joint rescaling of the x, t coordinates and 
the p coordinate; such a rescaling leads to a constant shift in /io,o- Thus, we can think 
of /io,o determining the scale of the cascading theory; note that this is independent of 
the temperature, since in the parametrization (4.2) all IR effects (including the effects 
of the temperature) are suppressed by powers of p. 

We would like to match (4.3) with the asymptotic solution (3.1) used above. We 
require that as p ^ (and correspondingly a; ^ 0) all the corresponding warp factors 
in the metric should agree to leading order, i.e., 



lim p-'hp?kp)-''' _ , kpf'kp) _ . 

{p,x}^o {1 - xy{2x - x^)-y^h{x)-^/^ ' {p,x}^oh{xy/^f2{x) 

{p,x}^Q h{xy/'^ f3{x) ' {p,x}^Qg{x) ' {p,x}^qK{x) 



(4.4) 
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This matching uniquely identifies: 

1 - 1 

^ = ^(^Ip^ + higher orders , 9q = go , /io,o Oq = ^0,0 + g^^^o ^^(y ) • (4-5) 

4.2 Expectation values in the black hole background 

In order to proceed, we would like to compute the expectation values of various op- 
erators in the cascading theories in terms of our parameters; in particular we want 
to compute the expectation value of the stress-energy tensor and of scalar operators 
which have dimension four when P — > 0. The expectation values of these operators 
were evaluated in [30] using the coefficients appearing in the expansion in powers of p 
of the functions appearing in (4.2), up to order p^. In order to recycle those results we 
need to translate our boundary expansion of the previous section to the one of [30], 
namely to write the leading terms of the expansion of [30] in terms of our parameters 
{/io,o, oo, 5*0, 02,0}, as we did for the zeroth order terms in (4.5) ^. 

When we do the matching we have some freedom, since in the parametrization (4.2) 
there is a freedom of performing diffeomorphisms of p depending on higher powers of 
p, that only affect the higher order terms in the expansion. Of course this freedom 
does not affect the eventual expectation values. Wc will fix this freedom by making an 
explicit choice for function of p to order 0{p^), given by : 

^ = \pX - p'al (^^ao + ^a^.o) + 0{p'') . (4.6) 

We can now identify all the terms in the expansion of [30] using our expansion of the 
previous section. Translating (4.5) to the notation of [30], we find at order p° 

1 1 a2 

po = ^0 , i^o = 4/io,oao - ^P'^Qo - ^P'^Qo In(^) , 

eg) =77,, = diag(-l,l,l,l) . 

All the coefficients in [30] at order p^ vanish, while the independent parameters ap- 
pearing at order are given by (again in the notation of [30]) 



a^'''^ = , a(^'^) = , a(^'^) = , gJ^^ = al , Gg? = , 

^Notc that the parameters {03 0,04,0} only show up at orders 0{p^) and 0{p^), respectively, so 
they do not affect the expectation values of these operators. 
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The one-point function of the stress energy tensor is given by [30] 

87rG5(T.,) = -icg) + 2GSf ^ + ^Gjf - a(^'°)) , (4.9) 

where G5 is the five dimensional Newton's constant obtained after we do the dimen- 
sional reduction on T^'^ (see [30]). In the normalizations that we are using, G5 is 
related to the ratio P/M (using the careful analysis of [35]) by = 87r"^P*^/81M^. 
Using (4.9) we obtain that the energy density £ and the pressure V are given in terms 
of our parameters by 

^ ,^ , 1 /3 2 6 
£ = {Ttt) = + -aoa2,o 

^^^^ ^ ^ , (4.10) 

^ = ^^--^ = 8^ " 7"°'''°. 

Since we do not have any chemical potentials, the free energy density J-' is 

The expectation values of the remaining two scalar operators which have dimension 4 
when P — are (using their normalization defined in [30]) 

24aoa2.o 



^^^"^ 7P^go 

n„) = 2^ + =^ (l + ln(^ 



10. . / .2 X (4.12) 

72, 

' 7go \' ' 2 

Note that in general curved backgrounds there was an ambiguity in some of the one- 
point correlation functions computed in [30], but there is no such ambiguity when the 
asymptotic four dimensional metric is flat (as in our case). 

4.3 The basic thermodynamic relation 

The equations above tell us, using the asymptotic values of the fields, that 

2 

sT = £-J^ = £ + V^-^ (4.13) 

in the cascading background, where s is the entropy density. On the other hand, we 
can also compute the entropy density and temperature directly at the horizon in terms 
of the horizon parameters {/iq, flo> fl'o > ^1} ■ 



4G5 ' 47r/iS6g V 4 + 2a1 ' ^ > 
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At first sight the previous two equations seem to give a non-trivial relation between 
some of our UV parameters and some of the IR parameters related to the expansion near 
the horizon. However, it turns out that sT is a renormalization group flow invariant 
in supergravity black brane geometries without a chemical potential [36, 37], so this 
relation is trivially satisfied in any solution of our equations of motion. 
To simplify notations we rewrite the metric (2.4) as 

2 

dsl^ = -cl de + cl {dxl + dxl + dx^'^+cl {dxf + cl {el)+cl (e^^ + e^J , (4.15) 



a=l 



where Cj = Cj(x) can be identified by comparing (2.4) and (4.15). Now, from the 
relation between components of the Ricci tensor 



R.7 



(4.16) 



we have a constraint^ 



4 4 
C2C4C5 

C3 



— 1 = constant 

C2, 



(4.17) 



Evaluating the left-hand side of (4.17) near the horizon, using the standard relations 
between the area of the horizon and the entropy and between the surface gravity of 
the horizon and the temperature, we have 



lim 



4 4 
C2C4C5 

C3 



SttG. sT . 



(4.18) 



On the other hand, evaluating the left-hand side of (4.17) near the boundary and using 
the asymptotic solution (3.1) we find 



lim 



4 4 
C3 



lim 

x^o+ - x^)G. 



1/2 - ^"0 ■ 



(4.19) 



Thus, our equation (4.13) follows in a straightforward way from the equations of motion. 



5 The numerical procedure 

5.1 Reducing the number of peirameters 

Before we begin the numerical solution of the equations, we can use the symmetries 
of the problem to get rid of some of our parameters. First, as mentioned above, the 

^Equation (4.17) can also be directly derived from (2.8)-(2.12). 
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parameters Q and ag are not scale invariant, but only the combination ho^al, which 
is a function of the temperature divided by the dynamical scale. We will choose as 
our parameter which is related to the dimensionless temperature the combination kg 
defined by 

P^goks = 4ho,oal - ^P^go ■ (5.1) 

Equation (4.5) now tells us that [kg — ln{al/2)/2\ is independent of the temperature 
(it depends only on the dynamical scale of the cascading theory). Thus, we can choose 
to define the scale A of this theory by a relation of the form 

Using the expressions for the high temperature entropy density of the theory computed 
in [11, 10], we sec that at high temperatures kg — (1/2) ln(T'*/A'*), with corrections 
scaling as ln(ln(T/A)). We will use kg instead of the temperature as our basic dimen- 
sionless parameter, and use (5.2) to translate between kg and T/A. 

Having understood this relation, we can now use the scaling symmetry (4.1) to 
set ao = 1, or, equivalcntly, use the fact that the solution to the equations of motion 
depends on the 7 UV parameters {go, ho^o, do, ci2,o, 92,0,(^3,0,(14,0} that we used in our 
expansion only through the six invariant combinations 

{, _ 0-2,0 ^ _ (l'3,0 ^ _ ^4,0 
90 ,ris , 0^2,0 = , 03,0 = , ^4,0 = , 92,0 

Oo Oo Oo 

Recall also that we are solving the theory in the supergravity approximation, which 
includes only the leading order terms both in the gg expansion and in the curvature {a') 
expansion. When we neglect gg corrections, the action (and the equations of motion 
we wrote) does not depend separately on and g but only on the combination P^g. 
We can thus set g^ — 1, and recall that whenever we have a factor of we really 
mean P^go. Furthermore, when we neglect a' corrections, the action is multiplied 
by a constant when we rescale the ten dimensional metric by a constant factor (and 
rescale the p-forms accordingly), so that the equations of motion are left invariant; this 
transformation acts on our variables as 

h^X-^h, /2,3^AV2,3, K^X^K, g^g, (5.4) 

and it changes P by P — > AP. We can use this transformation to relate the solutions 
for different values of P (as long as we are in the supergravity approximation). Thus, 
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we will perform the numerical analysis for P = 1, and we can use (5.4) to obtain the 
solutions for any other value of P. 

As a test of our numerics, we can check if it reproduces the solution at high temper- 
atures which can be computed pcrturbatively (as was done at leading order in [10], and 
at higher orders in the appendix). This computation implies that the correct solution 
should obey at large ks (A.35)-(A.38): 

/ln(2) 1021\ 1 /167809 (ln(2))2 781 , , u]\ 1 /, 

= + 1800 j ^ + (tosooo - T - 1200 ^"(2) + ' j ^ + ^ ) ' 

92,0 -{-l + l ln(2)) 1 + Q ln(2) - l(ln(2))^ + C?l) 1 + O (kf) , 

(5.5) 

where the values of the various constants appear in the appendix, and we used the 
high-temperature relation between the value of K at the horizon, which we denote by 
K^, and our dimensionless parameter kg'- 

= p2^o [ks + \ ln(2) + O {k:')^ . (5.6) 
5.2 Our numerical method 

As described above, for a given value of the temperature (or of kg) we have four 
parameters controlling the behavior of our solutions near the UV. What we need to 
do is to find for which value of these four parameters the solution is regular near 
the horizon, and this will determine the correct vacuum expectation values for this 
value of the temperature. The most naive way to proceed would be to go over all 
possible values of these parameters, use these values to determine the solution near 
the boundary, integrate the equations of motion up to a; = 1, and sec if the solution 
there is regular or not. Unfortunately, we cannot integrate the equations analytically 
but only numerically, and when we integrate the equations near the horizon, numerical 
errors always generate modes that blow up at the horizon, so we cannot really obtain 
solutions that are regular at the horizon in this way. 

One alternative might be to perform the integration in the opposite direction - start 
from a general solution near the horizon, integrate the equations to the boundary, and 
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see for which values of our near-horizon parameters we find a regular solution at the 
boundary (with the correct KT asymptotics) . However, this suffers from the same 
problem, that numerical errors generate modes that grow near the boundary. 

Thus, we are led to a procedure where we integrate the equations both from the 
boundary and from the horizon towards the middle of the interval x = 0.5, and attempt 
to match a solution that we get by integrating from the boundary with a solution that 
we get by integrating from the horizon. After setting qq = 1, for a given value of ks, the 
UV behaviour (3.1) is determined by 4 parameters (related to operator VEVs) (5.3) 

{^2,0 , ^3,0 , ^4,0 , 5'2,o} • (5.7) 
The IR behaviour (3.19) is determined by 6 horizon parameters 

{hlala1X:k^:9o}- (5-8) 

Matching a UV solution and an IR solution to (2.8)-(2.12) at a; = 0.5 implies 10 con- 
straints (5 for matching the values of the functions, and 5 for matching their derivative). 
Notice that we have precisely the same number of constraints as necessary to uniquely 
determine all the UV and IR parameters ((5.7) and (5.8)) for a given vahic of kg. 

Since both the boundary a; = and the horizon x = 1 arc singular points of the 
differential equations (2.8)-(2.12), we integrate the differential equations (2.8)-(2.12) 
from X = 0.01 (for the boundary integration) and from y = 0.01 (for the horizon in- 
tegration). In the former case the initial conditions are specified by the asymptotic 
expansion (3.1) which we developed to order x^/^ (inclusive); in the case of the hori- 
zon integration the initial conditions are specified by the asymptotic expansion (3.19) 
to order y^^ (inclusive). The coefficients of these asymptotic expansions generalize 
the results presented in section 3, and are available from the authors upon request. 
The mismatch between the boundary and the horizon integrations is encoded in the 
'mismatch vector' Vmismatch, defined by 



^mismatch — (^b ~ ^/i; K + f2,b — f2,h, /2,6 + /2,/i) /s.b ~ f3,h, fsfi + /3,/i) 



K, - Kh, K', + Ki g, - gh, g', + g'^^ 



(5.9) 



x=y=0.5 



where the subscripts ^ or ^ correspond to functions {h, /2, /s, K, g} integrated from the 
horizon or boundary, respectively, and the prime denotes derivatives with respect to x 
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Figure 1: Mismatch of hb{x) (for x < 0.5) and hh{y = 1 — x) (for x > 0.5) for different 
values of the parameters, for kg = 0.4. The sohd (blue) curves correspond to "correct" 
values of the parameters (5.7) and (5.8), with | It/mismate/il | ~ 9 x 10^^. The dotted 
(green) curves corresponding to all values of parameters 10% larger than the correct 
ones, produce | I^Viismate/il I ~ 3 x 10^^ The dashed (red) curves correspond to all values 
of parameters 20% smaller than the correct ones, giving | liVnismate/il I ~ 8 x 10~^. 

or y. The UV parameters (5.7) and the IR parameters (5.8) are tuned to ensure that 



An illustration of the integration as a function of the UV and IR parameters is presented 
in figure 1. 

We performed the numerical integration using Wolfram Mathematica® 6 with 40 
digit precision, to ensure sensitivity to the irrelevant operator parameters cis^o and 0,4,0- 

5.3 The numerical results 

We present the numerical results for the UV (5.7) and the IR (5.8) parameters as a 
function of kg in two regimes^°: 

-"^^The IR parameters are presented only for small values of fcs- Additional data are available from 
the authors upon request. 
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(5.10) 
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Figure 2: Values of the UV parameters 0-2,0 and 03,0 as a function of kg (blue points). 
The dashed/dotted (red/green) curves represent the perturbative 0{k~^) /Oik'"^) 
asymptotics of the parameters, given by (5.5). 

^4,0 92,0 

" ' 0.00 - 




Figure 3: Values of the UV parameters 04,0 and g2fl as a function of ks (blue points). 
The dashed/dotted (red/green) curves represent the perturbative 0{k~^) /Oik'"^) 
asymptotics of the parameters, given by (5.5). 

■ for large values of ks-, where we can check our numerical results against the pertur- 
bative analytic predictions (5.5); 

■ for an interval of small values of kg that includes the first order transition point 
to a confined thermal cascading background with broken chiral symmetry (as we will 
discuss in the next section). 

5.3.1 Large values of kg 

Figures 2 and 3 present the dependence of the UV parameters {02,0, 03,0, 04,0? fl'2,0} on 
kg G (4.29,24.0), with a step of Afc^ = 0.01 (blue points). In this regime the typical 
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Figure 4: The values of the UV parameters 02,0 and 03,0 as a function of ks 
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Figure 5: The values of the UV parameters 04,0 and (72,0 as a function of kg. 



norm of the mismatch vector (5.9) | l-ymismatc/il | ~ 10^^° or less. The dashed and dotted 
(red and green) curves represent the perturbative 0{k~^) and 0{kj^) asymptotics 
(5.5), respectively. Note that the dotted (green) curves track our numerical data quite 
well in this regime^^. 



5.3.2 Small values of kg 

Figures 4-8 present the dependence of the UV {a2,0) ^^3,0; ^^4,0; ^72,0} and the IR {a^, 
b^, g^, 4, h^} parameters on fc, e (0.25, 0.48) with a step of Afc, = 0.01 (Afc, = 0.005 
near the transition) (blue points). In this regime the typical norm of the mismatch 
vector (5.9) | ItViismatchl I ~ 10~^ or less. A highly non-trivial check on our numerics is 
the consistency of the holographic flow invariant sT. The latter can be computed in 

^^In appendix A we confirm using the perturbative high-temperature expansion that the cascading 
geometry thermodynamics satisfies the first law of thermodynamics up to order 0{k~^). 
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Figure 6: The values of the IR parameters Oq and function of 



Figure 7: The values of the IR parameters and as a function of kg 



the IR using (4.14), or in the UV using (4.13). We find that 

sT\iR 



sT\ 



- 1 



UV 



10' 



(5.11) 



or less, which provides an independent check on the accuracy of matching the IR and 
UV solutions (5.9). 

The special (red) point in figure 4 denotes a critical value of k^ = kcnucau for which 
the corresponding value of the parameter 0,2,0 = 0.2 flik critical) = ^ leads to a vanishing 
of the free energy density (4.11). We find kcriucai by performing a linear fit of the first 
5 numerical points: 

kcriucai = 0.25712(1) . (5.12) 
Our available numerical data shows that the free energy density (4.11) 



3al 



7 



28nG, \^'''~12 



(5.13) 
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Figure 8: The values of the IR parameters and as a function of kg. 

is negative when ks > kcriucai and is positive when kg < kcriucai, so we find that at 
temperatures lower than kg = kcriUcai the black hole solutions are not thermodynami- 
cally preferred over the gas of particles in the background of [9]. Thus, if we assume 
that these are the only two possible configurations, kcriUcai gives a critical temperature 
corresponding to a first order confinement transition, with chiral symmetry breaking, 
in the gravitational dual to the cascading gauge theory. Examination of the infrared 
parameters in figures 6-8 shows that the geometry at this transition is non-singular 
(as expected for a first order transition), and can be made arbitrarily weakly curved 
for large values of P, justifying the validity of the supergravity approximation. This 
observation is the main result of our paper. 

6 The physical results obtained from our numerical solutions 

In this final section we translate the results of the previous section into physical quan- 
tities in the theory as a function of the temperature. We present all the results as a 
function of T/A, where the temperature T is given by (4.14), and the scale A enters 
through the temperature dependence of kg. Recall that the numerical results presented 
in section 5 were obtained when setting P = qq = 1 and oq = 1. It is easy to restore 
the correct powers of P using the scaling symmetry (5.4), and to then put a factor of 
Qi) = 9o together with every factor of P^. In order to relax the Oq = 1 condition, all 
the dimensionful quantities must be computed in units of (see (5.2)) 
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In particular, from (4.14) we have 



A Airh^h'^ y 4 + 2a'l ' ^ ' ' 

enabling us to translate the dependence on kg into a dependence on the temperature. 

Equations (4.10) and (4.11) imply that the simplest expressions arise for the free 
energy density and the energy density divided by sT, which are given by 

Equations (5.2) and (4.14) allow us to compute the entropy density divided by the 
temperature cubed, which is a measure of the number of degrees of freedom in the 
theory : 

47rG'5 g _ 327r^gr _ 327r^ f _ ( ^ hi^h^H)^ - 9^)^ 
P^gl ~ 81M4^2 - 81M4^2 ^4 [j J ~ \^47r/i^6g y + 2a? 

' (6.4) 

Notice that at high temperatures we can use the perturbative expression (A. 42) of 
Appendix A and (4.13) to determine 

i_ = £l = ^ = ^ f 1 + o (^)) ^M*gl ^ (6.5) 

T3 T4 47rG'5T4 647rG'5 V \ J J 128 ^° ^A^ ^ ^ 

Finally, we can evaluate the vacuum expectation values of the two dimension 4 scalar 
operators (4.12) : 

(Oko) _24 e''^^ . 

(^Po) _.k. 92,, ^ 12(l-ln(2)+2A;,) ^ ^ ' ' 



Figure 9 presents In(^) as a function of kg at low and high temperatures. This is 
useful to determine the temperature dependence of the various UV and IR parameters 
presented in figures 2-8. Notice that the high temperature dependence of kg is in a good 
agreement with the high temperature asymptotic analysis of appendix A. Indeed, a 
straight line fit of the points in the plot on the right determines the slope to be 0.46(3), 
while the /c^ — > oo slope is expected to be | (A.40). 
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Figure 9: The relation between kg and the temperature T. 
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Figure 10: The free energy density JF, divided by sT, as a function of ^. On the left 
we plot temperatures at and slightly above the deconfinement transition, and on the 
right much higher temperatures. 



Figure 10 presents ^ as a function of ^ at low and high temperatures. Using a 
straight line fit of the first 6 points in the (left) free energy density plot, we determine 
the deconfinement and chiral symmetry restoration temperature to be 

=0.614111(3), (6.7) 

/ critical 

by requiring that the free energy density vanishes at T = Tcriticai- Notice that there are 
noticeably large deviations from scale invariant thermodynamics even for rather large 



temperatures. Indeed, for ^ ~ 10, the deviation of ^ from the conformal result 



conformal 



1 
1 



(6.^ 



is about 12%. 
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Figure 11: The energy density £, divided by sT, as a function of ^. 
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Figure 12: The temperature dependence of the effective number of degrees of freedom 
in the strongly coupled cascading gauge theory, as defined by N^fj oc s/T^. 

Figure 11 presents ^ as a function of ^ at low and high temperatures. Here, the 
deviation at high temperature from the conformal result 



S 



3 
4 



(6.9) 

conformal 

is three times less than the corresponding deviation in the free energy density (or 
pressure). Such a suppression is easy to understand once we notice from (6.3) that 



sT 



1 



sT 



sf 



I + -5] X — 
3 / sT 



(6.10) 

conformal \ / "'-^ conformal 

where 5 = Ya2,o- Note that the lattice results for QCD also imply that the energy 
density of the QCD plasma near the deconfinement transition is much more similar to 
that of scale-invariant thermodynamics than the QCD pressure [38]. 

Figure 12 presents the temperature dependence of the effective number of degrees of 
freedom of the strongly coupled cascading gauge theory at low and high temperatures. 



30 



4.0 



3.6 



T 

A 



0.620 0.625 0.630 0.635 0.640 



P 



I A4 



T 

A 



0.620 0.625 0.630 0.635 0.640 



Figure 13: The temperature dependence of the vacuum expectation values of the di- 
mension 4 operators {Okq) and {Opg). The operators are normahzed such that they 
are invariant under the scahng transformation (5.4). 



as defined by N^fj oc s/T^. One way to characterize the phase transition temperature 
is by the effective number of degrees of freedom (in the deconfined phase) at this 
temperature. Using the straight hne fit of the first 6 points in the (left) effective number 
of degrees of freedom plot, we find that at the deconfinement and chiral symmetry 
restoration temperature (6.7) 

327r^ i 



= 3.4291(5). (6.11) 

'-^critical 

Figure 13 presents the temperature dependence of the vacuum expectation values 
of the dimension 4 operators (Okq) and (Op^) (see (6.6)) at low temperatures. At 
high temperatures we can use the perturbative expressions (A. 42), (A. 35) and (A. 38) 
of appendix A to determine 



9o 



A4 
A4 



oc 



T 

In(^) 



In^ 



T 

'A' 



(6.12) 
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A Appendix : Perturbative evaluation of the solutions at high 
temperature 

In this appendix we analyze the high temperature thermodynamics of the cascading 
gauge theory perturbatively in P^qq/K^, where 



is the five-form fiux evaluated at the horizon, and kg is defined by (5.1). 

The purpose of this perturbative analysis is twofold: first, we would like to test our 
asymptotic identification of the cascading geometry parameters (4.5), of the temper- 
ature and of the dynamical scale A, against the first law of thermodynamics; second, 
we would like to obtain analytic predictions for the high-temperature values of the UV 
(5.7) and the IR parameters (5.8) perturbatively in P'^go/K-^, in order to benchmark 
our general numerical analysis. We will test the first law of thermodynamics to order 
O ^^^^ j inclusive, and evaluate the UV parameters (5.7) to order O (j-j^^ inclusive. 

The appendix is organized as follows. In subsection A.l we derive the perturbative 
equations of motion to order O (^-j^^ inclusive. In subsection A. 2 we present the 
near horizon and the near boundary expansions of the solutions, outline our numerical 
method for the computation of the UV/IR parameters of the perturbative solutions, 
and collect numerical expressions for some of these parameters. In subsection A. 3 we 
present perturbative expressions for the thermodynamics of the deconfined cascading 
gauge theory, and verify the first law of thermodynamics. 




(A.l) 
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A.l Perturbative equations of motion 

As discussed above, without loss of generality we can set go — go — 1. We use the 
following parametrization for the solution in perturbation theory in ^ : 



n=l 



f2(x) = So + ao 5^ I (^^^ (^-2^2n(x) + r]2n(x) + ^X2n{x)^ | + O (^^^ 



(A.3) 



fsix) =ao + ao ^ | (^^^ (^-2^2nix) + r]2nix) - ^\2nix)^ I + ao C 



p8X 



KtJ ' 
(A.4) 



= + E I (1^) «:2n(x)| + X. O (^1^) , (A.5) 

^(^) = 1 + E{(£)"c2.(a:)}+0(|^) . (A.6) 



The advantage of this parametrization is that the equations for {C2n,V2n, X2n,C2n} 
decouple, once the (decoupled) equation for K2n is solved, at each order {n — 1, 2, 3) in 
perturbation theory. We find (for n = 1, 2, 3) 



0^4n + ^ + , (A.7) 

""^2"+ ,r(x-l)(x-2) 3''2''^" + ^« ' ^^-^^ 
= + ^ - ^.^^ - 2.^ + , (A.10) 

= CI + ^ + 24 4n + ' 

3/ -L 

where the source terms {J^^ , jf-"^ , jf'^ , jf^"^ , J^f } are functionals of the lower 
order solutions: /t2m, ^2m^ V2m, X2m, (2111, with m < n. Explicit expressions for the 
source term functionals are available from the authors upon request. 
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The perturbative solutions to (A.7)-(A.ll) must be regular at the horizon, and 
must have the appropriate KT asymptotics (3.1) near the boundary. 

The leading order (n = 1) solution to (A.7)-(A.ll) was obtained in [10]: 

K2^--ln{2x-x^), (A.12) 
6 = -^ln(2x-a;^), (A.14) 

2 1 1 1 

C2 = ^ - 2dilog(l - x) + -dilog(2 - x) - - ln(x) ln(l - x) . (A.15) 

There is no simple expression for A2 — it is straightforward to write an appropriate 
solution using the Green's function for (A. 10), but this explicit expression is not useful. 

Similarly, although the higher order n = 2, 3 solutions to (A.7)-(A.ll) could be pre- 
sented in quadratures, these expressions are not useful. Rather, we identify the higher 
order solutions by specifying their asymptotic expressions near the horizon and near 
the boundary, along with the numerical values for the relevant integration constants. 

A. 2 UV/IR asymptotics of the perturbative solutions 

The UV/IR parameters of the higher order perturbative solutions are found by solving 
the differential equations (A.7)-(A.ll) numerically from the boundary {x — 0), and 
requiring the proper boundary conditions at the horizon, which are 

hm = lim rj^^ = hm ^'^^ = hm X'^^ = hm (L = , n = 2, 3 . (A. 16) 

X— S— >l_ >l_ X—*l_ X— 

To begin, we present the asymptotics of A2. As x ^ 0+ we find 

(A.17) 

where A3 ' is related to the condensate of the dimension 6 operator at order O ^ j . 
Nonsingularity of the X2{x) solution to (A. 10) at the horizon, together with 

hm A^ = , 

determines 

Af] = -0.872358024(9) . (A. 18) 
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Near the horizon, as y = 1 — a; ^ 0+, we find 

( 3181 5913 \ 8 ^ ln^ 
V 49152 16384 ^'V ^ ' ' 

where A20 can be determined numerically to be 

A^o = 0.16806(9). (A.20) 



A. 2.1 Order n — 2 asymptotics 



■ As X ^ 0+ we find 

+ lAr + , 

(A.21) 

(A.22) 

^4 = - ln(a;) + ( --k^^ + — - - ln(2) + - ln(a;)^ x - — A^'l x'/^ 
^ 36 ^ ^ V 6 ^ 144 24 ^ ^ 24 ^ 7 225 ^ 

+ f _ _ ^ J . . _ ^ 1^2/2) + (^il _ 1 in(2)^ ln(a;) (A.23) 

V 43200 36 ^ 80 ^ ^ 144 ^ ^ \720 12 ^ ' J ^ ' ^ ' 

_ ^ ln2(x)Yr2 - ^Af x^/^ + e» , 
144 ^ 7 3150 3 ^ V ; ' 

A4 = f-- + -4^1 + - ln(2)') X + aW + f — + + - ln(2) - ^ Inx") x' 

(A.24) 
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C4 = (cf'+(-|| + «?' + iM2))lnW):^ 



X 



(A.25) 

where the new UV parameters {^2^^ ^4^^, Ag^', (^'^} are determined numerically from the 
horizon boundary condition (A. 16): 



and 



K 



[4] 



^^3 



[4] 



0.73675974(3) , 



(A.26) 



numeric 



0.0053421556(6) , 



numeric 



= -1.1156300100(2) , 



(A.27) 



C: 



[4] 



numeric 



numeric 



0.622262593(4) 



With (A.26) and (A.27) there are no additional UV parameters to tune in order to 
enforce the horizon boundary condition (A. 16) for ^4. We find 



^'^{x = 0.99999) 



10" 



(A.28) 



numeric 



which is of the same order of magnitude as the error in (A. 16) for all other functions. 
■ The asymptotic expressions of the n — 2 solutions near the horizon y = 1 — x — > 0+ 
take the form 

«4 =4,0 + 4,2?/" + 4,4y^ + 0{y^) , 

^4 =^to + ^i2y' + ^'Ly' + 0{y') , (A.29) 

A4 =xio + A-^y" + Ad + (^(?/') , 
C4 ^Cjo + C4'2y' + + o{y') . 



To verify the first law of thermodynamics to order O ^ j we will need the numerical 
expressions only for {fi;^ 0, C4 O' ^4 2}- ^^^^ 



4,0 



= 0.62226(3) , 



(A.30) 



numeric 



36 



and 



?4,0 



?4,2 



numeric 



= -0.079819(3) , 

(A.31) 

= 0.019198(8) . 



numeric 



A. 2. 2 Order n = 3 asymptotics 

To verify the first law of thermodynamics to order O (^'^^^ we will need the asymptotic 
expression for only. We find 



1^6 = I + 



I - - I ln(2)) ln(.)) . + (1 + ^A^) .^/^ 



(A.32) 



/281821 2587 mi 509 , 31 x 1 m 13 m 1 J41 
+ (2I6OOO - 3600'^' - 720 ^^^2) + 288 ^^^^^ + 2^' + T^^' + 4^'' 

+ i ^-3^0 + S^^' - I ^^^2)) ^^^^) - lo 

f n53_ p] _ £ p] [4] 1 [4] 1 [2] A 5/2 ^ .3) 

V10500 ^ 25 ^ ^ 50 ^ 25 ^ ^ 7 ^ ^ 

as X — > 0+, with 

4^1 = -0.62226259(3) (A.33) 
determined from the horizon boundary condition (A. 16). 

A. 2. 3 Perturbative expressions for UV parameters (3.1) 

Finally, we collect perturbative expressions for the various independent UV parameters 
fl2,0i ^3,0) a,4,0i fi'2,0 as defined by (3.1). Because of the scaling symmetry (4.1) it is 
convenient to quote these parameters relative to ag- Also, in the next subsection we 
show that the first law of thermodynamics requires that 

4'' = §-^ln(2). (A.34) 

Note that (A.34) agrees with (A. 26) up to an error of order 10~^°. In the following 
expressions for the UV parameters in (3.1) we use the analytic expression (A.34). 
We find 

^"^o(^V (A.35) 



ao 12 
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(A.37) 

/I 1 \ P2 r 1 P"^ / \ 

9.o = (-5 + 5ln(2)) _ + cl"^ + o(^). (A.38) 

A. 3 Perturbative thermodynamics of the non-extremal cascading geome- 
try 

One of the interesting properties of the deconfined cascading geometry is the tempera- 
ture dependence of the five- form flux evaluated at the horizon K^, — Ki,{T). The need 
for such dependence was first pointed out in [10]; it stems from the fact that when 
studying the thermodynamics of non-conformal gauge theories (such as the cascading 
gauge theory) one must keep the intrinsic scale of the cascading gauge theory fixed^^, 
rather than keeping fixed the five-form at the horizon. 

The fact that K^, is temperature dependent introduces additional temperature de- 
pendence into the thermodynamic potentials (the free energy density (4.11), the 
energy density S (4.10), and the entropy density s (4.14)) via the UV parameters oq 
and 02,0, both of which depend on K^,. As a result, the first law of thermodynamics 

dJT = -s dT (A.39) 

would not be valid, unless the temperature dependence of is properly determined 
and taken into account. One possible approach is to use the first law of thermodynamics 
(A.39) as a way to determine fC^(T). Such an approach was proposed and implemented 
in [30] to leading order in O , where it was found that validity of (A.39) requires 
that 

,KAT)_2P^^^(^^ (A.40) 



dT T \K^^ 
(this was also proposed in [10], based on the requirement of keeping the glueball spec- 
trum scale of the cascading gauge theory fixed). Equation (A.40) was also shown to 



^This fact was not clearly taken into account in previous numerical studies [13]. 
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be required for the consistency of the hydrodynamics of the cascading gauge theory 
plasma in [18]. 

The main observation of this paper is that one can rigorously determine the tem- 
perature dependence of iiT* without referring to the first law of thermodynamics. In 
the bulk of the paper this was implicitly done in our solutions. In the context of the 
perturbative high-temperature expansion, we can obtain such an identification pertur- 
batively by expanding the exact matching condition (4.5), enforcing the fixed scale 
of the cascading gauge theory, perturbatively in We will demonstrate here that 
this identification is consistent with the first law of thermodynamics. This provides a 
non-trivial consistency check on our solutions. 

In the rest of this subsection we present explicit expressions for function 
of the temperature T to order O ^^j- can then use (A. 35) to compute the 

thermodynamic potentials of the cascading black hole geometry. We present explicit 
perturbative expressions for '^^^^^ , and verify the first law of thermodynamics to order 

"(«)■ 

A. 3.1 Cascading black hole thermodynamics to order O 

Explicitly evaluating the temperature of the black hole to order O (^-^^ we find 

«.^^(i+l^.c(nv 



4 V V^.' 
Using (A. 12)- (A. 15), we further determine 



4 \ 2K^ 
The matching condition (4.5) then determines 

0(g^4V„-^. + P'^(iK^^)-i). (A.43) 

Assuming that — K^{T) and differentiating (A.43), we find (A. 40). 

A. 3. 2 Cascading black hole thermodynamics to order O (^^^ 

We can evaluate the black hole temperature to order O j by requiring that the 
Euclidean continuation of the metric (2.4) does not have a conical singularity asy ^ 0+. 
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We find^^ 

2 f-aoY'j. ^P' f in' + 8)' 1 , 2 , , 2 

^^\k:) \ Kl \ 1920 + 30 ^^2'°^ ^^'^ + ^^^'^ 3 ''^'^ 

*'(fi)} 

(A.44) 

Solving for ao from (A.44), and reading off {ao, 02,0} in (3-1) from (A.2)-(A.6), we find 

1^2 2.. f P"^ / 13 1 , , , TT^ TT^ 1 , , ,2 4 

+ 2eto-4et,)+o(|^)}, 



(A.45) 



7^2 -202 f. 



Additionally we find (see (5.1)) 



P^k. ^ Ah„,^l - ip^ = if.{l-MH)^ + o(gj|, ( A.47) 

We are now ready to verify the first law of thermodynamics. The matching condition 
(4.5), to order O {jr^^ gives 

a (I) = 4.„,„ - i.. + (1 .„(ll^) - i) + ^ . (A.48) 

which results in the following ordinary differential equation for = '■ 



dT T { \K^^ , 

Now, given (A.45) and (A.46) we can evaluate the energy density £ and the pressure 
V. The first law of thermodynamics (A. 39) leads to 

Consistency of (A. 49) and (A. 50) makes a prediction 

4''-§-^ln(2). (A.51) 



i^We used from (A.19) X!^ ^ = ("i + |^2,o)- 
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As a highly non-trivial check on our numerical analysis of the perturbative expansion, 
note that (A. 51) agrees with (A. 26) to within a factor of order 10~^°. 
With (A. 51) we can also evaluate the speed of sound squared 



2 dV 1 4 PV 10 PV 



A. 3. 3 First law of thermodynamics to order O 

For the temperature dependence of to order O we can again find two expres- 

sions — one involving both the UV parameters and the IR parameters, and the other 
one parameter-independent. These parallel the expressions (A. 50) and (A. 49) : 

^ =^ + ^ + ^ (- 480^ - 30^ + 9 ^^^2) +180-15 ^^^^) + 



(A.53) 

dK, 2P2 f P2 1 P4 /P6\1 



dT T { 2 Ky 

where we used (A. 34). Consistency of (A.53) and (A.54) leads to a prediction 



1 



2 X 10-^ . (A.56) 



-io"'-^"'+^^"^')+^-^ (Ay%^<o+<o-8et2+cr+2Kr' = 1. (A.55) 

We can estimate the error in our solutions by comparing the two sides of (A.55). Using 
the explicit expressions for the perturbative UV parameters (A. 27) and the perturbative 
IR parameters (A. 30) and (A. 31), we find 

LHS 
RHS 

Finally, the speed of sound can be expressed either in terms of the UV parameters 
or the IR parameters 

"3 9 K,^ 27 Kl + ^ 81 9 9^^ ) Kl^^\Kt) ' ^ ^ 

2 1 4 P2 10 P^ / 1 ^ 2 2 4 , 11 8 .2 

^^^3-9 K: + 27KF+i-T080" "135" +81^"(')- 45 - 135 (^^'°) 
32 ^ IQ 1^ _ 32 1^ \ P6 / p8 \ 

(A.58) 

Consistency of (A. 57) and (A.58) is guaranteed by (A.55). 
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